Ground State Properties of Simple Elements from GW Calculations 
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A novel self-consistent implementation of Hedin's GW perturbation theory is introduced. This 
finite-temperature method uses Hartree-Fock wave functions to represent Green's function. GW 
equations are solved with full potential linear augmented plane wave (FLAPW) method at each 
iteration of a self-consistent cycle. With our approach we are able to calculate total energy as a 
function of the lattice parameter. Ground state properties calculated for Na, Al, and Si compare 
well with experimental data. 

PACS numbers: 71.10.Ca, 71.15.Nc. 71.20.-b 



Density-functional theory (DFT)[1] in its local density 
approximation (LDA) or generalized gradient approxima- 
tion (GGA) is a widely used method to calculate ground 
state properties of solids. But this theory is not always 
good in situations where electronic correlations are essen- 
tial. The problem is then arising that existing variants 
of DFT cannot be improved systematically (or at least 
it is very difficult to accomplish that). It is therefore 
very desirable to have a method capable to deal with 
ground state properties and at the same time allowing 
their systematical improvement. This flexibility is obvi- 
ously present in methods based on diagrammatic expan- 
sion and one of them, Hedin's GW approach, is becoming 
computationally accessible during last years. However, 
while the GW method has become standard for studying 
excitation spectra, the possibility of using it for calcu- 
lating total energies and other ground state properties is 
not well established. This question has been addressed 
only in a few works. Holm and von Barth [2| have per- 
formed self-consistent GW calculations for the homoge- 
neous electron gas (HEG) and concluded that in spite of a 
bad electronic structure obtained in their self-consistent 
(SC) calculations the total energy was quite close to the 
result of a quantum Monte Carlo Q study. Garcia- 
Gonzales and Godby [H also have obtained good total 
energies of HEG (including spin-polarized case) in their 
SC GW calculations. Stan, Dahlen and van Leeuwen 
Q have applied SC GW to calculate total energies of 
atoms and molecules. They conclude that GW calcu- 
lations should be done self-consistently in order to ob- 
tain physically meaningful and unambiguous energy dif- 
ferences. Miyake, Aryasetiawan, Kino, and Terakura @ 
have studied ground state properties of sodium and alu- 
minum using Galitskii-Migdal formula with a model 
spectral function. Their equilibrium volumes and bulk 
modulus appeared to be slightly overestimated as com- 
pared to the experiment, but some improvement over 
LDA results was reported. Miyake et al. [8| have applied 
the total energy formula due to Luttinger and Ward [9( to 
the calculation of equilibrium lattice constants in Na and 
Si. Their one-shot type results appeared to be very close 
to the experimental data. However, to our knowledge 



fully self-consistent GW calculations of the total energy 
and the ground state properties for real solids have not 
yet been carried out. 

To address this question we have implemented a vari- 
ant of GW method which allows us to calculate total en- 
ergies. The key ingredients of our implementation are the 
following. First, we use full potential linear augmented 
plane waves method (FLAPW) to find the solutions 
of the Hartree-Fock (HF) equations which serve as a ba- 
sis for the expansion of one-electron Green's function. 
Second, we have found it vital to use Mazubara's time 
r-mesh to calculate correlated part of the self-energy. 
It appears that it is very difficult to obtain comparable 
accuracy in total energy using self-energies calculated 
in frequency domain. Third, our calculations are self- 
consistent. 

In our imlcmentation there are two formulas for the 
total energy. They differ in how to find its exchange- 
correlation part. The first implementation uses the 
Galitskii-Migdal formula, i.e. convolution of the Green 
function G\i\ (k; t) and the self-energy £™ A , (k; — r) 
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where a is a spin index; k denotes points in the Brillouin 
zone, A, A' are band indexes. Second formula follows from 
the equation of motion for the one-particle Green func- 
tion: 
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where Eki n is a kinetic energy of the electrons; Vwj^k) 
denote matrix elements of the Hartree potential; /i is a 
chemical potential; (3 = l/T, and Q is the number of 
electrons. 

We check the convergence of the total energy with re- 
spect to the number of points on r-mesh by comparing 
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FIG. 1: Total energy of HEG with respect to r s obtained in 
HF and GW approximations. Comparison is made with QMC 
results. |3| 



TABLE I: Parameters of calculations 





Na 


Al 


Si 


k-mesh 


8x8x8 


8x8x8 


5x5x5 


bands 


23 


15 


40 


N T 


80 


48 


40 


Nlapw 


23 


38 


73 


N PB 


135 


190 


513 


Semicore 


2s2p 


2s2p 





the result from Eq. ([2|) with the result from Eq. (p} . We 
find that the formula ([1]) gives much more stable results 
than Eq.©. 

As a test of our newly developed code we have applied 
it to calculate the total energy for the homogeneous elec- 
tron gas. As we have already mentioned, this was al- 
ready done previously (@, [3), but in those calculations 
zero temperature formalism was utilized. In Fig. [T] the 
result from our finite temperature (300K) approach is 
shown in comparison with the accurate Quantum Monte 
Carlo data obtained by Ceperly and Alder. Q We have 
also calculated total energy of HEG using the HF ap- 
proach. As it is seen, SCF GW total energies are in very 
good agreement with QMC simulations and show great 
improvement compared with the HF results. 

All our calculations for real solids have been performed 
for temperature T=2000K. For our LDA and GGA stud- 
ies we used th e exchange-correlation functionals from 
Refs. 



11] 



[121 ] respectively. We have listed most im- 
portant calculational parameters for all studied systems 
in Table [J 

They include the k-point sampling in the Brillouin 



zone, the size of LAPW basis, Nlapw, the size of prod- 
uct basis, Npb, the number of bands used for the rep- 
resentation of correlated part of the Green function, and 
the number of points on Mazubara's time mesh. Abso- 
lute convergence of the total energy in our calculation 
is a bit worse than the one typical for LDA calculation 
(lmRy/atom), but it should also be noted that the con- 
vergence of DOS is reached faster than the convergence 
of the total energy. Furthermore, there is some depen- 
dence of our predicted equilibrium volumes with respect 
to some specific parameters. The most critical one is the 
number of k-points. The influence of other parameters 
mostly results in a rigid shift of the total energy with no 
much change in the position of a minimum. 

As an example, Fig. [5] presents our calculated total 
energy versus lattice parameter for Na while Table |TT] 
contains the calculated ground state properties for all 
studied elements: Na, Al, and Si. Here we compare the 
results of the LDA, GGA, HF, and GW approximations. 
All these calculations are self-consistent. As it is seen 
there is the same tendency for all materials with respect 
to the predicted lattice parameters: the LDA gives too 
small values, while the HF approach predicts them to be 
too big (especially for metals). The GGA result is bet- 
ter, though for Al this approach leads to a bit expanded 
volume. We should stress that the GW result is very ac- 
curate (especially for aluminum and silicon). For sodium, 
both the GW and the GGA results are quite similar. We 
should also note that for a semiconductor (Si) the HF 
approach is quite competitive in accuracy with the GW 
but it gives a little expanded volume which is not good 
if one thinks of additional expansion connected to the 
lattice vibrations. 

Bulk moduli do not show such a consistent improve- 
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FIG. 2: Total energy of Na as a function of lattice parameter. 
Arrow indicates the experimental lattice parameter. [131] 
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TABLE II: Equilibrium lattice parameter ao(a.u.) and bulk 
modulus Bo(GPa) of Al, Na, and Si compared to experimental 
data [H, [l| 



Na. 



Al 



Si 





ao 


Bo 


ao 


Bo 


ao 


Bo 


LDA 


7.53 


87 


7.585 


82.5 


10.01 


91.2 


GGA 


7.83 


71.5 


7.72 


73.3 


10.14 


93.5 


HF 


8.62 


50 


7.82 


82.5 


10.34 


97.0 


GW 


7.82 


63.0 


7.64 


86.5 


10.17 


100.7 


Exp 


7.96 


68.1 


7.65 


72.16 


10.26 


99 



merit when calculated within the GW approach. Only for 
silicon we have found a remarkable agreement with ex- 
periment. We think that the reason is that our results are 
not well convergent with respect to the k-points sampling 
used but it is too expensive computationally to study its 
full convergence. 

As a byproduct of our study we have investigated the 
influence of the self-consistency on the excitation spec- 
tra. Despite it was already noticed long ago that 
the GW approximation can be quite useful for calcu- 
lating one-electron spectra, there is still some ambigu- 
ity with respect to the effect of self consistency. Usu- 
ally, GW calculations are exceedingly demanding and 
only one-shot variant is used to calculate band widths 
or band gaps. Calculated on top of LDA or GGA 
such one-shot quasiparticle band structures appeared to 
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FIG. 3: DOS of Al calculated in GGA, HF and GW approx- 
imations. Chemical potential is placed at zero energy. The 
arrow shows the experimental position of the valence band 
bottom. 



be in much better agreement with experimental data 
when compared to the LDA ones 0, Gil, [H, H, 
Later it was shown that self consistency within GW 
gives too large band gaps for semiconductors and insu- 
lators 17 ,[3, [2l|. Shishkin, Marsman and Kresse 21 1 



_ |21|. Shishkin, Marsman and Kresse 
have argued recently that the inclusion of vertex cor- 
rections is required to obtain accurate band gaps in the 
framework of SC GW. However this question of self- 
consistency still remains as a few other works aimed 
at self-consistent GW calculations used some simplifica- 
tions. Schone and Eguiluz 1171 utilized pseudopotential 
approach. Ku and Equiluz [22| have applied all-electron 
method but with 'diagonal' approximation to Green's 
function. The calculations by Zein, Savrasov, and Kotliar 
[23| were performed with atomic sphere approximation 
(ASA). Bruneval, Vast, and Reining [3] have combined 
the self-consistent screened exchange plus Coulomb hole 
(COHSEX) with one-shot GW. Shishkin and Kresse Q 
have updated self-consistently only the eigenvalues in the 
Green's function. Kotani, van Schilfgaarde, and Faleev 
[20| | applied a strategy of finding an effective Hamilto- 
nian so that further perturbative contributions would be 
as small as possible. Their approach has some justifica- 
tions though it is not true answer for the self-consistent 
GW. 

To address the issue of an influence of the self- 
consistency on the calculated electronic structure we per- 
form SC GW calculations and compare the results with 
the results from the one-shot GW calculations based on 
DFT self-consistent calculations. As example of calcu- 
lated DOS for metals, we present our calculated elec- 
tronic structure of Al in FigJ3] We are especially con- 
cerned with the valence band width for this material. 
First observation (which may be trivial one) is that HF 
approximation gives too big band width. Second obser- 
vation is also quite known for simple metals: both LDA 
and GGA calculations give slightly expanded band width 
as compared to experimental data. As it is seen, our 
non SC GW band width for Al is much closer to the 
LDA result and to experimental data than SC GW re- 
sult which is about 10% too big as compared to the ex- 
perimental works by Lyo and Plummer[l3j and by Livins 
and Schnatterly [25j. Since Aluminum is a free-electron 
metal, we think that the above result is consistent with 
the results by Holm and Barth 0], and by Shirley [26| 
obtained for the homogeneous electron gas, where the 
calculated by SC GW band width was found to be too 
big as compared to the non self-consistent one. Based on 
the work by Shirley we expect that higher order vertex 
corrections will bring the calculated band width in closer 
agreement with experiment. 

In Table HIIl we have collected the calculated fundamen- 
tal gaps for silicon obtained in earlier works along with 
our results. Our band gap from non SC GW calculation 
is quite close to the results of others. However our band 
gap from SC GW appears to be a bit wider than it is 
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TABLE III: Band gap (energies in eV) for Si. The values in round brackets were obtained with SC in eigen values only. 
Experimental band gap is 1.17ev. 



Ref.[17] Rcf.[22] Rcf.[23] Ref.[20] Ref.[19] Rcf.[21] Present work 



LDA 0.53 0.52 0.46 0.51 0.46 

GGA 0.54 

HF 6.27 

GldaWlda 1.34 0.85 0.86 0.98 1.14 0.86 

GhfWhf 2.69 

GcohsexWcohsex 1.56 

GWlda 1.28(1.20) 

QPscGW 1.25(1.14) 1.47 1.41 

GW 1.91 1.03 1.10 1.55 



usually obtained. We would stress however that neither 
of previous calculations are approximation free and that 
partially SC results obtained in the works [3 M, El 
also show a trend in increasing the band gap. In this 
respect it would be very interesting to know what effect 
will bring the vertex correction on the band gap in Si. 
This question we hope to answer in the nearest future. 

In summary, we have presented a self-consistent real- 
ization of the GW method and its performance for the 
evaluation of the total energy and ground-state prop- 
erties. For the materials we studied (Na, Al, and Si) 
the GW approach leads to consistently better descrip- 
tion of equilibrium volume than LDA, GGA or HF ap- 
proaches do. Our results for the electronic structure of 
the above mentioned materials are in qualitative agree- 
ment with earlier works, but in general we should con- 
clude that some deterioration in calculated band widths 
and gaps is seen when we are trying to use all-electron, 
full-potential method in connection with " classical" self- 
consistent GW approach. In this respect our conclusion 
is the same as one by Holm and Barth 0] : self-consistent 
GW approach can produce good total energies but it is 
not so accurate for the one-electron spectra. 
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